Search for: All records

A new version of the US National Science Foundation National Center forAtmospheric Research (NSF NCAR) thermosphere-ionosphere-electrodynamicsgeneral circulation model (TIEGCM) has been developed and released. Thispaper describes the changes and improvements of the new version 3.0since its last major release (2.0) in 2016. These include: 1) increasingthe model resolution in both the horizontal and vertical dimensions, aswell as the ionospheric dynamo solver; 2) upward extension of the modelupper boundary to enable more accurate simulations of the topsideionosphere and neutral density in the lower exosphere; 3) improvedparameterization for thermal electron heating rate; 4) resolvingtransport of minor species N(2D); 5) treating helium as a major species;6) parameterization for additional physical processes, such as SAPS andelectrojet turbulent heating; 7) including parallel ion drag in theneutral momentum equation; 8) nudging of prognostic fields near thelower boundary from external data; 9) modification to the NO reactionrate and auroral heating rate; 10) outputs of diagnostic analysis termsof the equations; 11) new functionalities enabling model simulations ofcertain recurrent phenomena, such as solar flares and eclipse. Wepresent examples of the model validation during a moderate storm andcompare simulation results by turning on/off new functionalities todemonstrate the related new model capabilities. Furthermore, the modelis upgraded to comply with the new computer software environment at NSFNCAR for easy installation and run setup and with new visualizationtools. Finally, the model limitations and future development plans arediscussed. 

Abstract This paper investigates the midlatitude ionospheric disturbances over the American/Atlantic longitude sector during an intense geomagnetic storm on 23 April 2023. The study utilized a combination of ground‐based observations (Global Navigation Satellite System total electron content and ionosonde) along with measurements from multiple satellite missions (GOLD, Swarm, Defense Meteorological Satellite Program, and TIMED/GUVI) to analyze storm‐time electrodynamics and neutral dynamics. We found that the storm main phase was characterized by distinct midlatitude ionospheric density gradient structures as follows: (a) In the European‐Atlantic longitude sector, a significant midlatitude bubble‐like ionospheric super‐depletion structure (BLISS) was observed after sunset. This BLISS appeared as a low‐density channel extending poleward/westward and reached ∼40° geomagnetic latitude, corresponding to an APEX height of ∼5,000 km. (b) Coincident with the BLISS, a dynamic storm‐enhanced density plume rapidly formed and decayed at local afternoon in the North American sector, with the plume intensity being doubled and halved in just a few hours. (c) The simultaneous occurrence of these strong yet opposite midlatitude gradient structures could be mainly attributed to common key drivers of prompt penetration electric fields and subauroral polarization stream electric fields. This shed light on the important role of storm‐time electrodynamic processes in shaping global ionospheric disturbances. 

Abstract The Global‐scale Observations of Limb and Disk (GOLD) and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instruments were used to investigate the thermospheric composition and temperature responses to the geomagnetic storm on 23–24 April, 2023. Global‐scale Observations of Limb and Disk observed a faster recovery of thermospheric column density ratio of O to N₂(ΣO/N₂) in the southern hemisphere (SH) after the storm ended at 12 Universal time (UT) on 24 April. After 12 UT on 25 April, ΣO/N₂had mostly recovered in both hemispheres. Global‐scale Observations of Limb and Disk also observed an increase of middle thermospheric temperature (140–200 km) (Tdisk) on 24 April with a maximum of 340 K. Within 4–6 hr of the storm ending on 24 April, Tdisk enhancement persisted between 30°N and 60°N, 100°W and 30°W, while Tdisk lower than pre‐storm quiet day (17 April) was observed between 45°W and 15°W, 40°S and 50°N. Tdisk recovered between 100°W and 45°W, 30°N and 55°S. On 25 April, Tdisk was lower than on 17 April across the entire GOLD Field‐of‐Regard (FOR) by ∼50–110 K. Additionally, solar irradiance decreased by 15%–20% from 17 to 25 April, indicating that the lower Tdisk on 25 April resulted from both storm and solar irradiance variations. Latitudinal variations of Tdisk and the SABER observed Nitric Oxide (NO) cooling rate revealed that NO cooling is crucial for the lower Tdisk in the northern hemisphere (NH) mid‐high latitudes on 25 April. These results provide direct evidence of decreased thermospheric temperature during storm recovery phase than pre‐storm quiet times. 

Abstract ICON observations were used to investigate local time (LT) and latitudinal variations of thermospheric meridional winds in the middle‐high thermosphere (160–300 km) during quiet times in 2020 June and December. At middle‐low latitudes (10°S–40°N), meridional winds were predominantly equatorward in the summer hemisphere while mostly poleward in the winter hemisphere. The meridional winds showed that the diurnal variation was dominant between ∼20°N and ∼40°N, but the semi‐diurnal variation played a leading role at lower latitudes (below ∼20°N) during solstice months. Thermosphere‐Ionosphere Electrodynamics General Circulation Model reproduced the ICON observed meridional wind variations qualitatively. A model diagnostic analysis shows that the pressure gradient force dominated the semi‐diurnal variation of the winds, while the Coriolis force played a leading role in the diurnal variation in June. In December, LT variations of meridional winds were primarily driven by pressure gradient and ion drag forces. During both months, the vertical viscosity was important, tending to balance the effects of pressure gradients. Additionally, semi‐diurnal variations of low‐latitude meridional winds in June were more affected by upward propagating tides than those in December. 

Abstract This work investigates mid‐ and low‐latitude ionospheric disturbances over the American sector during a moderate but geo‐effective geomagnetic storm on 13–14 March 2022 (π‐Day storm), using ground‐based Global Navigation Satellite System total electron content data, ionosonde observations, and space‐borne measurements from the Global‐scale Observations of Limb and Disk (GOLD), Swarm, the Defense Meteorological Satellite Program (DMSP), and the Ionospheric Connection Explorer (ICON) satellites. Our results show that this modest but geo‐effective storm created a number of large ionospheric disturbances, especially the dynamic multi‐scale electron density gradient features in the storm main phase as follows: (a) The low‐latitude equatorial ionization anomaly (EIA) exhibited a dramatic storm‐time deformation and reformation, where the EIA crests evolved into a bright equatorial band for 1–2 hr and then quickly separated back into the typical double‐crest structure with a broad crest width and deep equatorial trough. (b) Strong equatorial plasma bubbles (EPBs) occurred with an abnormally high latitude/altitude extension, reaching the geomagnetic latitude of ∼30°, corresponding to an Apex height of 2,600 km above the dip equator. (c) The midlatitude ionosphere experienced a conspicuous storm‐enhanced density (SED) plume structure associated with the subauroral polarization stream (SAPS). This SED/SAPS feature showed an unusual temporal variation that intensified and diminished twice. These distinct mid‐ and low‐latitude ionospheric disturbances could be attributed to the storm‐time electrodynamic effect of electric field perturbation, along with contributions from neutral dynamics and thermospheric composition change. 

Abstract We report a new ionosphere phenomenon: Equatorial ionization anomaly (EIA) discontinuity (EIAD), based on OI 135.6 nm radiance observations from the Global Observations of Limb and Disk (GOLD), ground‐based total electron content maps and in‐situ ion density data from Constellation Observing System for Meteorology, Ionosphere, and Climate‐2. The EIAD occurs when the OI radiance of the EIA crest has a local minimum, at a fixed UT, with the radiance in the local longitude region being weaker than that on the east and west sides. In the GOLD field‐of‐view, EIAD follows the seasonal variations of EIA. EIAD appears more often over the Atlantic Ocean and Africa than over South America. It occurs more in the southern crest during the December solstice, and more in the northern crest during both equinoxes. EIAD can occur under both quiet and disturbed times. 

Abstract The geomagnetic storm on February 3, 2022 caused the loss of 38 Starlink satellites of Space‐X. The Global‐scale Observations of the Limb and Disk (GOLD) observations and Multi‐Scale Atmosphere Geospace Environment (MAGE) model simulations are utilized to investigate the thermospheric composition responses to the Space‐X storm. The percentage difference of the GOLD observed thermospheric O and N₂column density ratio (∑O/N₂) between the storm time (February 3, Day‐of‐Year [DOY] 34) and quiet time (DOY 32) shows a depletion region in the local noon sector mid‐high latitudes in the southern hemisphere, which corresponds to the east side of GOLD field‐of‐view (FOV). This is different from the classic theory of thermospheric composition disturbance during geomagnetic storms, under which the ∑O/N₂depletion is usually generated at local midnight and high latitudes, and thus, appear on the west side of GOLD FOV. MAGE simulations reproduce the observations qualitatively and indicate that the ∑O/N₂depletion is formed due to strong upwelling in the local morning caused by strong Joule heating. Interestingly, enhanced equatorward winds appear near local midnight, but also in the local morning sector, which transports ∑O/N₂depletion equatorward. The depletion corotates toward the local afternoon and is observed in the GOLD FOV. The equatorward winds in the local morning are due to the ion‐neutral coupling under the conditions of a dominant positive interplanetary magnetic field east‐west component (B_y) during the storm. 

Abstract Following the 2022 Tonga Volcano eruption, dramatic suppression and deformation of the equatorial ionization anomaly (EIA) crests occurred in the American sector ∼14,000 km away from the epicenter. The EIA crests variations and associated ionosphere‐thermosphere disturbances were investigated using Global Navigation Satellite System total electron content data, Global‐scale Observations of the Limb and Disk ultraviolet images, Ionospheric Connection Explorer wind data, and ionosonde observations. The main results are as follows: (a) Following the eastward passage of expected eruption‐induced atmospheric disturbances, daytime EIA crests, especially the southern one, showed severe suppression of more than 10 TEC Unit and collapsed equatorward over 10° latitudes, forming a single band of enhanced density near the geomagnetic equator around 14–17 UT, (b) Evening EIA crests experienced a drastic deformation around 22 UT, forming a unique X‐pattern in a limited longitudinal area between 20 and 40°W. (c) Thermospheric horizontal winds, especially the zonal winds, showed long‐lasting quasi‐periodic fluctuations between ±200 m/s for 7–8 hr after the passage of volcano‐induced Lamb waves. The EIA suppression and X‐pattern merging was consistent with a westward equatorial zonal dynamo electric field induced by the strong zonal wind oscillation with a westward reversal. 

Abstract The extension of the neutral sodium (Na) layer into the thermosphere (up to 170 km) has recently been observed at low and high latitudes using a Na lidar. However, the geophysical mechanisms and implications of its formation are currently unknown. In this study, we conduct an advanced two‐dimensional numerical simulation of the Na and Na⁺variations in theEandFregions at low latitudes. The numerical simulations are used to investigate the contributions of the electromagnetic force, neutral wind, diffusion, and gravity. The simulations lead to three major findings. First, Na⁺in the subtropical region of the geomagnetic equator acts as the major reservoir of the neutral sodium, and its distribution during nighttime is mostly below 200 km due to the combined effect of the vertical component of thedrift and Coulomb‐induced drift. Second, we find that the fountain effect has little influence on the behavior of Na in the nighttime. Third, the probable explanation for the frequent generation of the thermospheric sodium layer during spring equinox at Cerro Pachón, Chile is attributed to the large vertical neutral transport generated by large vertical wind perturbations of unknown origin, with a magnitude exceeding 10 m/s that is closely associated with the semidiurnal tide. 

Abstract In order to understand the characteristics of long‐lasting “C‐type” structure in the Sodium (Na) lidargram, six cases from different observational locations have been analyzed. The Na lidargram, collected from low‐, middle‐, and high‐latitude sites, show long lifetime of the C‐type structures which is believed to be the manifestation of Kelvin‐Helmholtz (KH) billows in the Mesosphere and Lower Thermosphere (MLT) region. In order to explore the characteristics of the long‐lasting C‐type structures, the altitude profile of square of Brunt‐Väisälä frequency in the MLT region has been derived using the temperature profile collected from the Na lidar instruments and the SABER instrument onboard TIMED satellite. It is found to be positive in the C‐type structure region for all the six cases which indicates that the regions are convectively stable. Simultaneous wind measurements, which allowed us to calculate the Richardson numbers and Reynolds numbers for three cases, suggest that the regions where the C‐type structure appeared were dynamically stable and nonturbulent. This paper brings out a hypothesis wherein the low temperature can increase the magnitude of the Prandtl number and convectively stable atmospheric region can cause the magnitude of Reynolds number to decrease. As a consequence, the remnant of previously generated KH billows in nearly “frozen‐in” condition can be advected through this conducive region to a different location by the background wind where they can sustain for a long time without much deformation. These long‐lived KH billows in the MLT region will eventually manifest the long‐lasting C‐type structures in the Na lidargram.